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1. Introduction 

Finite Size Scaling (FSS) provides a powerful tool which is used to extrapolate information to 
infinite volume. It can also be used to study the critical behavior and calculate critical exponents. 
When the system is close to a 2nd order phase transition, some long distance properties can be 
determined by the global features of the system, such as the symmetries and the dimensionality. 
Universality enables us to classify seemingly different systems into certain classes. Svetitsky and 
Yaffe [P pointed out that, for d + 1 dimensional SU (N) pure gauge model, "if any portion of the 
boundary is second-order, then the critical behavior will be described by some fixed point of d- 
dimensional, Z(N) invariant spin systems." For a pure gauge SU (2) theory in 3+1 dimension, we 
expect the universality class of the finite temperature transition to be the same as the 3 dimensional 
Ising model. Existing results on FSS for SU (2) [Q, pL ffl, pi [q] agree well with this expectation. 

In the following, we work on SU (2) gauge theory in 3+1 dimensions: N T x N„, where N r = 4, 
and N a = 2, 4, 6, 8, ■ • • 16, with periodic boundary conditions. We focus on direct estimations of the 
exponents v and CO. For N t = 4, the only direct estimate of v we are aware is 0.65(4) obtained 
in Ref. [g], which is compatible with the accurate value 0.6298(5) obtained in Ref. [0] for the 
3-dimensional Ising model. We would like: 1) to improve the accuracy of the estimate of Ref. [||] 
and 2) resolve the corrections due to the irrelevant directions. Part 2) is largely unexplored and 
a better understanding of these corrections could help us design methods to reduce these effects 
as done in Ref. [0] for the 3 dimensional Ising model. In the existing work for = 4 in Refs. 
[0, ||, Pi, a fixed j3 interval procedure was used. This means the j3 interval is fixed for different 
volumes. In these proceedings, we shrink the interval in order to reduce the nonlinear effects [^] 
and use a finer j8 resolution. 

2. Binder cumulants and FSS 

In the following, we define the 4th order Binder cumulant [Q], g$, as 

(p4\ j i N r 

Related definitions appear in the literature, such as B4 = jpiji- We assume that there is no external 
field and that the g4 depends on the scaling variables as 

g 4 = g 4 {u K Nc /v ,u l N a m ,...) (2.2) 

with 

u K = K + u^'K +... 
U\ = uf^ +U\^K+ ... 

where u K is the only relevant scaling variable, u\ is the first irrelevant scaling variable, and K is the 
reduced quantity K = (j3 — j8 c ) / j8 c . Expanding up to the first nonlinear corrections, we obtain 

^ 4 (j8,A^ CT ) =g4(PcH+fiKN l J v + f 2 K 2 N* /v + (c + c l KN l J v )N a <a + --- (2.3) 
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This expansion is accurate if \k\No V is small enough. In addition we would like the nonlinear 
effects of fi and c\ to be negligible. It is easy to estimate fi from numerical data for intermediate 
values of K ^ . However, c\ is more difficult to resolve from an already small effect and its effect 
will be ignored. In the following, we will work with values of | k\n]J v such that the effects of fi 
are within the numerical errors of #4 that we now proceed to discuss. 

3. Error Analysis 

For each volume, we worked on the region near j8 c (we discuss j8 c later). We generated 50,000 
Polyakov loops with 20 different seeds for each j3. Therefore, for each j8, we have 1,000,000 values. 
However, the data is correlated and we will need to remove the correlations. The autocorrelation 





Figure 1: These two graphs show that the autocorrelation time changes with the volume. For N a =A, the 
autocorrelation time is around 10. For ^=16, the autocorrelation time is around 200 illustrating the critical 
slowing down. 
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Figure 2: variation of P and P from the ensemble average of for 10 different seeds . 



times change with j8 and the volume. Fig. [T] shows two autocorrelation functions. In order to 
remove the autocorrelation among the data, we tried the bootstrap and jackknife methods. For 
the jackknife method, we skipped every T data in order to remove the autocorrelation, where T is 
the autocorrelation time. Comparisons with other methods to estimate the errors will be discussed 
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elsewhere. We should note here that one should not calculate the error of P 4 and P 2 first and 
then use the error propagation to get the error of g 4 . This is simply because P 4 and P 2 are highly 
correlated, see Fig. ||. The jackknife method can remove the correlation between P 4 and P 2 . 

4. Determination of the critical exponent v 

We now focus on the estimation of v. If we are reasonably close to j3 c , we can use the linear 
form: 

g 4 (P,N a ) ^g 4 (j8o-) + co^V ff a) + /i^y v (4.1) 

This expression contains 6 unknown parameters: gA(Pc,°°), cq, (0, f\, f5 c and 1/v and we will use 
a new strategy to attack this difficult problem [|l0|]. A first observation is that the dependence on v 
can be isolated from the other parameters by studying the linear dependence in j3. However, one 
should keep in mind that the slope of this linear fit depends implicitly on the choice of the center of 
the interval which should be as close as possible to f3 c .. In order to guarantee that nonlinear effects 
are under control, we will consider j3 intervals of the form 

1/3 - (p c ) app \ < 0.015 x (4/N a f^ (4.2) 

We will start with reasonable values for (P c )app and (l/v) app and then show that the effect of their 
variations is small. The factor 0.015 has been chosen following a procedure described in Ref. [^] 
and guarantees an approximate linear relation between g4 and j3 at fixed N a - 

g A ^a Na +b Na x j8 . (4.3) 

This is illustrated in the left side of Fig. g for {p c ) app = 2.299, (1 /v) app = 1.6 and N a = 10. We 
see that in the chosen interval the deviation from linearity seem mostly due to statistical errors 
rather than a systematic curvature that would be observed if we had chosen a broader interval. The 
volume dependence of the parameters of the linear fit are 

a Na ~ g 4 (p c , oo) + coN a m - f Y N l J v (4.4) 
b Na ^fiN l J v (4.5) 

Once we have for different N a , we can do a log-log fit to determine the inverse critical exponent 
1/v. One should note that the slope can be determined independently of /3 C or g4(j3 c ,°°). This is 
illustrated in the right side of Fig. g. From the log-log fit with N a > 6, we obtain 1/v = 1.56(4). 
We plan to model and explain the deviations from linearity at low N a . We need to address the 
dependence on (j3 c ) and (l/v) app . We used the central value f5 c = 2.2991 from Ref. [|]] as 
the critical value and changed the center of the interval (j5 c ) ap p, between 2.297 and 2.301 and 
(V V W between 1.4 and 1.8. We calculated 1/v for a set of 189 values of ((j6 c ) app , 0-/v) app ). 
The histogram is shown in Figure | left. We can see that the values appear in a rather large range 
between 1.468 and 1.631. The average of this set of results is 1/v = 1.570 with a = 0.027 which 
will be our preliminary result. At the conference, a larger value was found, but with low statistics at 
large N a . In the meantime, we collected more data at large N a . If we choose the center of the data 
in a smaller range, namely with Q3 C ) from 2.298 and 2.300 and the same range for (1/v) , we 
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Figure 3: Left: linear fit of #4 near /3 C for N a = 10, {fic)app = 2.299 and (1/v) = 1.6. We get a Nrj and 
bx a from the fit for different volumes. Right: determination of 1 / v from the log-log fit discussed in the text. 
bi ~ f\ x n]J v / /3 £ . By taking Log[\bN a |] vs. Log[N a ], we can get the slope which is just 1/v. We only did 
this fit with N a between 6 and 16. The Na = 2, and 4 points deviate from the linear approximation. 
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Figure 4: Left Fig.: {fi c ) apP changes from 2.297 to 2.301; {l/v) app changes from 1.4 to 1.8. l/v= 1.570; 
(7 = 0.027. Right Fig.: {fi c )app changes from 2.298 to 2.300; {l/v) app changes from 1.4 to 1.8. TJv = 1.571; 
(7 = 0.028. 

get a slightly different distribution shown on the right side of ||, but the average is essentially the 
same. It seems thus possible to average over (p c ) ap p and (1/v) in order to get a more accurate 



value of 1/v. For reference, l/v /ri is estimated as 1.5887(85) in Ref. [JXTJ] and 1.5878(12) in 
Ref. [Q] and our preliminary result here is consistent with these more accurate values. This is 
preliminary, we plan to go to larger and to determine j3 c independently using the method of Ref. 



5. Determination of the critical exponent 00 

Unless we determine j8 f and 1/v very precisely, it is very difficult to subtract the effects of the 



third term of Eq. ( p~T| ). If we can work at j8 c , this term is absent: 

g4(Pc,N a )= g4 (p c ,™)+coN a a 



(5.1) 



Consistently with the previous section and the rest of the literature, we assume the universal value 
g4{Pc°°) = 0.46575 as found in Ref. [f7|]. Log[\g4 — g4(j3 c ,°°)|] vs. Log[N a ] should be linear right 
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Figure 5: For J3 C =2.2991, the behavior is approximately linear: g4 ~ g4(/3 c ,°°) +coxN a m . 



at p c and nonlinear for all the other /3s. This is shown in Figure |5} At the same time, the slope 
00. The result we obtained from this analysis is oo = 2.030(36). This is very different from 



is 



ooismg = 0.812 [IT]1. It is possible that the coefficient of the N a a is very small and the exponent we 
extrapolated is a sub-subleading exponent. For a detail discussion of the subleading corrections, 
see Ref. [Q]. The most plausible explanation seems that this exponent is related to the irrelevant 
direction associated with the breaking of rotational symmetry [12] and which is close to 2. 



6. Conclusion and perspective 



By using methods discussed in Refs. [H, 10], we analyzed the 4-th order Binder cumulant 



of pure SU(2) lattice gauge theory and estimated the critical exponent v and oo. 1/v is in good 
agreement with the value for the universality class of the 3D Ising model, while the corrections to 
scaling seem dominated by anisotropic effects fll2[]. 
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